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ABSTRACT 

We calculate the Yukawa corrections of order 0{aew'm'^f^i^^/m'yy), C'(aeu,"7.j(.j-j/TO^) 
and 0{aew'mt{b)/''^w) ^'^ widths of the decays t2 ti + (/i°, in the Min- 

imal Supersymmetric Standard Model, and perform a detailed numerical analysis. 
We also compare the results with the ones presented in an earlier literature, where 
the 0{as) SUSY-QCD corrections to the same three decay processes have been cal- 
culated. Our numerical results show that for the decays t2 ti+h'^ , <2 ^ ti+H''\ 
the Yukawa corrections are significant in most of the parameter range, which can 
reach a few ten percent, and for the decay t2 ti+A'^ , the Yukawa corrections are 
relatively smaller, which are only a few percent. The numerical calculations also 
show that using the running quark masses and the running trilinear coupling At 
, which include the QCD, SUSY-QCD, SUSY-Electroweak effects and resume all 
high order (tan /3)-enhanced effects, can vastly improve the convergence of the per- 
turbation expansion. We also discuss the effects of the running of the higgsino mass 
parameter fi on the corrections, and find that they are significant, too, especially 
for large tan/3. 



PACS number(s): 14.80. Cp; 14.80.Ly; 12.38.Bx 



1 Introduction 



Incorporation of supersymmetry is one of the most attractive and promising pos- 
sibilities for new physics beyond the Standard Model (SM)|1|, 0], and the Minimal 
Supersymmetric Standard Model (MSSM) is a popular candidate for new physics in 
this way. In the MSSM there are many new particles. For example, every quark has 
two spin zero partners called squarks qi and qr, one for each chirality eigenstate, 
which mix to form the mass eigenstates qi and g2- For the third generation quarks, 
due to large Yukawa couplings, there may be large mass differences between the lighter 
mass eigenstate and the heavier one, which implies in general a very complex decay 
pattern of the heavier state. 

As we know, the next generation of colliders, such as the Large Hadron Collider 
(LHC), the upgraded Tevatron, e^e~ linear colliders, and collider will push the 

discovery reach for supersymmetric (SUSY) paticles with masses up to 2.5 TeV0, |^ 
and allow for precise measurement of the MSSM parameters. Thus a more accu- 
rate calculations of the deacy mechanisms beyond the tree level are necessary. The 
dominate decay modes of the heavier squarks are shown as below: 

All these squark decays have been extensively discussed at the tree-level[|, |, 0]. 
Up to now, one-loop QCD and supersymmetric (SUSY) QCD corrections to above 
decay channels have been calculated too P, ||, ^, while the Yukawa corrections and 
the full electroweak one-loop radiative corrections to the squark decays into quarks 
plus charginos/neutralinos were given in Ref.[|Tn] and Ref.||ll||, respectively. Also 



the Yukawa corrections to the squark decays into quarks plus gluinos were given in 
Refs. |p!2|, p!3| , and the Yukawa corrections to the heavier squark decays into lighter 
squarks plus vetor bosons were given in Ref. [Q. Recently, the Yukawa corrections 



to the bottom squark decays into lighter top squarks plus charged Higgs bosons has 



been presented in Ref . [|15] . So only the electroweak radiative corrections to the heavier 
top squark decays into hghter top squarks plus neutral Higgs bosons have not been 
calculated yet, including the Yukawa corrections to these processes. 

In this paper, we present the calculations of the Yukawa corrections of order 

0{aewfnt(b)/™w)' ^{(^ew^t{b)/''^w)y O {aewf^t{b) I "^w) widths of the heav- 

ier top squark decays into lighter top squarks plus neutral Higgs bosons, i.e. the de- 
cays ^2 ^1 + {h^ 1 1 A^) . These corrections are mainly induced by the Yukawa 
couplings from Higgs-quark-quark couplings, Higgs-squark-squark couplings, Higgs- 
Higgs-squark-squark couplings, chargino(neutralino)-quark-squark couplings, and squark- 
squark-squark-squark couplings. As shown in Ref. , the Higgs-Squark-Squark cou- 
plings receive large radiative corrections, which can make the perturbation calculation 
of the relevant Squark or Higgs boson decay widths quite unreliable in some region 
of the parameter space. When the correction term is negative, the corrected width 
can even become negative, which clearly makes no sense. In order to solve this prob- 
lem, we use the running quark masses and the running trilinear coupling At [|l6l , and 



vastly improve the convergence of the perturbation expansion. We also discuss the 
effects of the running of the higgsino mass parameter /i on the corrections, and find 
that they are significant, too, especially for large tan/3. 

2 Notation and tree-level result 

In order to make this paper self-contained, we first summarize our notation and 
present the relevant interaction Lagrangians of the MSSM and the tree-level decay 
rates for ts ^ ti + A^). 

The current eigenstates qi and qR are related to the mass eigenstates qi and q2 by 

\ (12 J \ 1R J \ -smOg cosOq J 

with < < TT by convention. Correspondingly, the mass eigenvalues m^^ and rrig^ 
(with < rrig^) are given by 

(1- ,ni)-^'^i^^'^'- z^:) 



with 



= M^ + ml + mlcos2p{IlL-e,sm^ew), (3) 
^Ir = ^{u,D} + T^l + ^z cos 2peq sin^ Ow, (4) 
Qq = Ag — fi{cot P,taii P} (5) 

for {up, down} type squarks. Here M| is the squark mass matrix. Mqjj I) and At^ 
are soft SUSY-breaking parameters and fi is the higgsino mass parameter . I^^^ and 



Cq are the third component of the weak isospin and the electric charge of the quark 
q, respectively. 

Defining Hk = {h°,H°,A°,G°,H^,G^) (fc=l,...,6), one can write the relevant 
lagrangian density in the {qi,q2) basis as following form (z,j=l,2; a and (3 are flavor 
indices) : 

Relevant = H^fiatPL + b'^Mq'^ + {Gt^Huqf'^qt + gq{4kPR + 

+ {Gf,),,HiH,q^*qr + h.c, (6) 

with 

(GfJ,,- = [R^G%{R^n, {I, k = 1, 6) (7) 



{Gt),, = [R^GURh% {k = l,...,6) (8) 

where and are the couplings in the (g^, qn) basis , and their explicit forms are 
shown in Appendix A. The notations a^, 6^ (k=l,...,6), and a^^,, bf^. (k=l,...,4), and 
^iky kit; (k=l,2) used in Eq.(|D are defined also in Appendix A. 

The tree- level amplitudes of the three decay processes, as shown in Fig. 1(a), are 
given by 



^ {At cos a + yU sin a) 



{At cos a + /i sin 



a 



-^/2mtht cos a + 



grriz sin(a+/3) 



G, 



tR 



{RThi (9) 



3 



for t2 



r(o) _ -iDt I -V^iTT'tht sin a - 9m,cos{a+i3) sin a - /i cos a) | ^ ^^-^^^ 

^(Asina-/icosa) -V2mtht sin a - ^^^^^^^^Ct, 



for t2 — ^ i^i-f^", and 

Mf = ^[i?*Y , + (11) 

^ 2mw \-Atcotp-fi y V ; J21 v ; 

for ^ Here h, = h = ^j^^,, C,, = - e..^, = e,4, 

svK = sin 6*^, and = cos6't„. I^j^ = \ and Ct = | for the top squark, = — | and 
Cft = — I for the bottom squark. The tree-level decay width is thus given by 

^ |M(o)pAV^K,m|,<o) 

where X{x,y,z) = {x — y — z)"^ — Ayz and s=(l,2,3) corresponds to the decay into 
respectively. 

3 Yukawa corrections 

The Feynman diagrams contributing to the Yukawa corrections to t2 — > tiH^ are 
shown in Figs.l(b)-(f) and Fig.2. We carried out the calculation in the t'Hooft- 
Feynman gauge and used the dimensional reduction, which preserves supersymmetry, 
for regularization of the ultraviolet divergences in the virtual loop corrections using 



the on- mass-shell renormalization scheme |]T^, in which the fine-structure constant 



Uew and physical masses are chosen to be the renormalized parameters, and finite 
parts of the counterterms are fixed by the renormalization conditions. The coupling 
constant g is related to the input parameters e, mw and mz via g'^ = e^/s^ and 

= 1 — m^/m|. As for the renormalization of the parameters in the Higgs sector 
and the squark sector, it will be described in detail below. 

The relevant renormalization constants are defined as 



mqo = rriq + Smq, m|o = m?+(5m?, 

AqO = Aq + 6Aq, = jj + Sfi, 

dqo = 0q + 66q, tan Pq = {1 + 6 Z fj) tau P , 

sinao = (1 + ^Za) sin a, 

q^ = il + 6Zfy/' + 6Z!^q,, 

H° = {1 + SZho^/^H^ + SZHOhoh^, 

/i° = (1 + 6ZhoY/^h^ + SZhOHoH^, 

Go = (1 + 5Zg-Y/'G- + 5ZghH-, 

AO = (1 + 5Zao)^/'A° (13) 

with q = t,b. Here we introduce the mixing of H~ and . 

Taking into account the Yukawa corrections, the renormahzed amphtude for t2 — ^ 
tiHg is given by 

where SM^''^ and SM^f^ are the vertex corrections and the counterterms, respectively. 
The calculations of the vertex corrections from Fig. 1(b)- 1(f) result in 

k=l j 

i ^ - - - 

fc=l ij 

+ Y.(^lk)jiiGl)2jBo{ml,mHO,mq^) 

k=l j 

-Y^Esin^?cos^,-(^?4^$i - hlRil^^2)iGlhBoimlo,m-^m-,J 

[3(sin^e,- + cos^0,-) -2sin2 0,-cos2^,-](G*)2ii?o(m|o,m,-^,m,-J 
-8 sin^ Of cos^ 6^(01) uBq (m^o , m^-^ , m^-^ ) 
+4 sin 9^cos9fCos29f{Gl)iiBo{m1[o, mi_^,mi^) 
-4 sin 6*^ cos Oi cos 26't-(G* ) 22^0 ("^^0 , m^-^ , m^-^ ) } 
+Fxs- (15) 



where F^s is the remains, which are given by 

ifif^/lt COS(a — (S — 1)^) ^ , , , r 7* ; J *, nr, O ^„ 

Fxs = - ' l^^l, [2(4al. K][(2m?^+P*.Pi^o)Cn 

+ i'2PhPH0 + m^o)Ci2 + 2m?^C2i + 2m^oC22 + ^Pi^PH^C23 + 8C24] 

+ [2(4,6*/ + 6*,a*/)m^o][(2m| +p?,p^^o)Cn + (2pt,p^^o + m^o)Ci2 
+ (m|C2i + m^oC22 + 2pi^pHoC23 + 46*24 + ("^1 + Pi^PHo)Co + m^Co] 
}iPii^PHO,m^o,mt,mt) 

+(2k,Pho + m^o)Ci2 + 2m|C2i + 2m^oC22 + ^Pi, PH2C23 + 8C24] 



+ [2(4^H + )^v-] [(2m4 + Pi,PHo)Cn + {2pi,PH0 + mU)a 
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+ (m|C2i + m^oC22 + '2pi^PHoC23 + 4C24 + (m| + Pi^PHo)Co + m^Co] 
}iPh,PHO,rn^-,mb,mb) (16) 
for s=(l,2), and 

i^x3 = -^^^ E{ [2(44; - V*]("^ioCi2 + Pi.PAoCn) 

+ [2(a2i^*ii* - ^2ia*ii*)"^x?][(2"^| + PhPA")Cn + C^Pi.PAO + m^„)Ci2 
+m?^C2i + m^oC22 + 2pi^pAoC23 + 4C24 + (p| + PhPAo)CQ - m^^Co] 
}(Pfi,PAo,m^o,mt,mi) 

+ [2(44* - 44*)^x-][(2< +PfiPAo)Cn + (2pi^p^o +mi„)Ci2 
+m|C2i + m^oC22 + 2pi^pAoC23 + 4C24 + {p\ + PtiPAo)Co - m^Co] 
}{Ph^PAO,m--,mb,mb) (17) 

for s=3. In above expressions i?o and are two- and three-point Feynman integrals pC 
respectively. For q = t, we have k = 1...4. For q = b, we have = 5, 6. 
The counterterms can be expressed as 

5Mf ) = tiG{)2i l^iSZi + 5Z2 + 5Zho) - 2 tan 29^59^ + — - - cos^ pSp] 
2 g 2my^ 



.gmtAtCosa AtSmt + rritSAt , , 2 

-I — cos 2u^[ h tan aoZa 

2mw sm p rritAt 



-I- — cos 2ei[— + 5Za\ 

2mw sm p n rrit 

+t{G{)uSZ^2 + tiGlh2SZ2i + tiGl)2iSZH0HO, (18) 



SM^'' = iiGl)2i [UsZi + 6Z2 + 5Zho) - 2 tan 26^66^ + — - - cos' p5/3] 

2 g 2my^ 



.gnitAtsma At5mt + mt5At , , 

-I- — - cos 2et[ + 6Za\ 

2mY/ sm p TTitAt 



.gumtcosa r^/^ , ^"^t 2 

I coszfe'f 1 tan i 

2mw sm p fi rrit 



+i{Gl)u5Z,2 + tiGl)225Z2i + tiG\)2iSZhOHO, (19) 



(5Mf = i{Gl)2i[-{5Z^ + 5Z2 + 5Zao) + 



2 



2 ' ' g 2m 



2 J 



cot/3[ 5Zp\-- [— + ]. (20) 



2raw rritAt 2mw ^ rat 

Here we consider only the counterterms from the Yukawa couphngs, and the exphcit 
expressions of some renormahzation constants calculated from the self-energy dia- 
grams in Fig. 2 are given in Appendix B. Other renormahzation constants are fixed 
as follows. 

For 6Zgh, using the approach discussed in the two-Higgs doublet model (2HDM) 
T8[| , we derived below its expression in the MSSM, where the version of the Higgs 



m 



potential is different from one of Ref. |[1^|. First, the one-loop renormalized two-point 
function is given by 

iTcHip^) = i{p^ - m]j-)SZHG + ip^SZcH - iTgh + iT^cnip'^), (21) 
where Tgh is the tadpole function, which is given by 

Tgh = tt—ITh, sin(« Th, cos(a - /?)]. (22) 

2mw 

Next, from the on-shell renormahzation condition, we obtained 



1 



8Zgh = ^^[Tgh - T.GH{m'H-)]. (23) 

The explicit expressions of T^gh and the tadpole counterterms Tu^ {k = 1, 2) are 
given in Appendix B. 



For the renormalization of the parameter /3, following the analysis of Ref . |T9[ , we 
fixed the renormalization constant by the requirement that the on-mass-shell H^fVr 



coupling remain the same form as in Eq. (3) of Ref. |T9| to all orders of perturbation 
theory. However, with introducing the mixing of H and G~ instead of H~ and W~ , 
the expression of 5Zp is then changed to 



2 



2 mi 2 mi — mi/ 2 



For the counterterm of squark mixing angle Oq, using the same renormalized 
scheme as Ref.[|TOi, one has 



56^ 



Re[Sf,K) + sL(m|)] 
2(m?^ -m|) 



(25) 



where the explicit expressions of the Sj^ functions arising from the self-energy dia- 
grams due to the Yukawa couplings are given in the Appendix B. 

For the renormalization of soft SUSY-breaking parameter Aq, we fixed its coun- 



terterm by keeping the tree-level relation of Aq, iriq. and 9q ||2T|, and get the expression 
as following: 



m| - m? ^ ^ ^ 5m„, sin 20,7 

"'(2 cos 29,69. - sin 29g — + 



2mn 



2mn 



q2> 



+{cot /3, tan/3}5/i + 5{cot /5, tan (3}^. 



(26) 



As for the parameter /x, there are several schemesfTll, ^ to fix its counterterm 



and here we use the on-shell renormalization scheme in Ref. [23 1, which gives 



6iJi = '^[m-+{6Uk2Vk2 + Uk26Vk2) + 6m^+Uk2Vk2], 



(27) 



k=l 



where (?7, V) are the two 2x2 matrices diagonalizing the chargino mass matrix, and 
their counterterms {6U, 6V) are given by 



6U=-{6Zr-6ZI)U, 



4 



6V = -{6Zl-6ZI)V. 



(28) 
(29) 
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The mass shifts Sm-+ and the off-diagonal wave function renormahzation constants 
SZFi{L) can be written as 

5m-+ = lRe[m-+(n^,(m^,) + n« (m|,)) + nff (m^,) + nff (mj,)], (30) 



■mt+ — mt+ ■' '^i ^] ■' 

Xi Xj 



+n^'"(mj+)m-+ + n^^"(m^+)m-+], (31) 



i5ZL)^, = {5Zn),, (L^R). (32) 

The exphcit expressions of the chargino self-energy matrices n^(^) and n'^'^(^) are 
given in Appendix B. 

Finally, the renormalized decay width is then given by 

r, = r(°) + 5T^:^ + (5ri^) (33) 

with 

^1/2 ( jjj^ 777-^ TDj^ 1 

^p(a) ^ """^ ReiMi^^^t^Mi^H (a = v,c). (34) 

4 Numerical results and conclusion 

We now present some numerical results for the Yukawa corrections to the decays 
h ti + (/i°,if°, A°). The SM input parameters in our calculations were taken to be 



aewimz) = 1/128.8, mw = 80.375GeV, mz = 91.1867GeV|2l, rrit = 175.6GeV, and 
rub = 4.25GeV. 

In order to improve the convergence of the perturbation expansion, using the 
method presented in Ref.|]l6|, we take into account the QCD and SUSY QCD run- 
ning quark masses mg{Q){mt{Q),rhb{Q)) and running trilinear coupling At in our 
calculation(the energy scale Q here is the mass of the heavier top squark i.e. wi^j)- ■'^^ 
the tree-level H^t2ti couplings, we use rht{Q) and At instead of the on-shell parame- 
ters. While in the calculation of the one-loop corrections, all parameters are on-shell 
except the Yukawa Couplings ht, hb taken as the running quark masses. 
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rhq{Q) are evaluated by the next-to-leading order formula p5|, p6[ 

rhb{Q) = U6{Q,mt)Urj{mt,mi,)mi,{mb), 

rnt{Q) = UG{Q,mt)mt{mt), (35) 

where we have assumed that there are no other colored particles with masses between 
scale Q and rrit, and rhbi^b) = 4.25GeV,mt(mj) = 175.6GeV[^. The evolution 
factor Uf is 

- J(/) - 8982 - 504/ + AOP 

33-2/' 3(33 -2/)2 ' 



where as{Q) is given by the solutions of the two-loop renormalization group equations[^ 
When Q = 400GeV, the running mass rhh{Q) ~ 2.5GeV. In addition, we also im- 
proved the perturbation calculations by the following replacement ^ 



* 1 + Amt{MsusYQCD)' 

•mbiQ) — — r 7TT r, 38 

1 + /\mi,[MsusY) 



where 



= - — {Bi{0,mg,mtj + Bi{0,mg,mtj - sin26'i( — )[5o(0, m^, m^-J - 5o(0,m 



2a. ^ /i? 



'^T^b = —Mg/i tan /5/(m^^,m^2,Mg) + — ^fiAt tan /3 1 {mi^^m^^, fi) 



167r2 



/iM2 tan /?[cos 6il{mi^ , M2, /i) + sin 9il{mi^ , M2, yu) 

+ ^ cos^ 9~J{mi^,M2, /i) + ^ sin^ ^j;/(m^^, M2, //)] (40) 



with 



C) = (,2_^2)(^2!,2)(,2_,2) (^^^^ + ^ + (41) 

The running trilinear couplings At can be obtained according to the procedure of 



the DR renormalization where the UV divergence parameter A = 2/e — 7 + log47r is 

10 



set to be zero |1I6|. First we compute the running stop masses m? (Q) = mj +5mj and 



the running mixing angle of the top squarks 6i{Q) = 0i + 56i, where the counterterms 
(5mf and 59; are given by 



561 = - ^ , ''''^ . (42) 

2 mi — mi 

Here the exphcit expressions of the Sj^ functions arising from the QCD self-energy 



diagrams are given in Ref||2g]. Then we can get the running parameter At from the 
formula 

rhtAt = {ml (Q) - m| (Q)) sin OiiQ) cos ^t(Q) + rht^i cot (3. (43) 



The two-loop leading- log relations 1 29] of the neutral Higgs boson masses and mix- 
ing angles in the MSSM were used. For mH+ the tree-level formula was used. Other 
MSSM parameters were determined as follows: 

(i) For the parameters Mi, M2 and fi in the chargino and neutralino matrix, we 
take M2 and fi as the input parameters, and then use the relation Mi = {5/3){g'^/g'^)M2 
O.5M2II, |30| to determine Mi. The gluino mass mg was related to M2 by mg = 
{as{mg)/a2)M2^. 

(ii) For the parameters ^ and At^b in squark mass matrices, we assumed 
Mq = 1.5Mjj = l.bMfj and At = At, to simplify the calculations, except for Figs. 10- 
11, where we assumed M^, = 1.12Mq and At = Af, in order to compare with the 
SUSY-QCD results in Ref.g. 

Some typical numerical results of the tree-level decay widths and the Yukawa 
corrections are given in Figs. 3- 12. 

Figs. 3 - 5 show the m^-^ dependence of the results of the three decay channels, 
respectively. Here we take m^o = ISOGeV, fi = M2 = 200GeV, and At = Af, = 
600GeV. The leading terms of the tree-level amplitudes Mj°)(s=l,2,3) are given by 

^ = (At cos a + u sin a) cos 201 , (44) 

2mw sm p 
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^ = !^!!^ (it sin a - u cos a) cos 201 , (45) 

2mvi/ sm p 



Mf = -^(Acot/5 + /i). (46) 

For mt^=100GeV, cos^^- ~ (-0.575, -0.574, -0.574) and cos a ~ (0.754, 0.953, 1.000) 
for tan/3 = 4, 10, and 30 , respectively, and for mt^=560GeV, cos^^f ~ (-0.323, -0.332, 
-0.334) and cos a ~ (0.737, 0.897, 0.992) for tan/3 ~ 4, 10, and 30, respectively. In 
the case of z = 2, the two terms in Eq.(^) have opposite signs, and their magnitudes 
get close with the increasing tan/3 and thus cancel to large extent for large tan/3. 
Therefore, the tree-level decay widths have the feature of ro(tan /3 = 4) > ro(tan (3 = 
10) > ro(tan/3 = 30) in most range of the parameter space, as shown in Fig.4(a). 
In the case of i = 1, the two terms in Eg . (^4]) have the same signs, there are not 



cancelling effects between them, so Fq is larger than the one of the case of z = 2 
for the same values of tan/3. In the case of z = 3, the amplitude contains a term 
propotional to cot /3, so Fo(tan/3 = 4) > Fo(tan/3 = 10) > Fo(tan/3 = 30). From 
Figs.3-5(b), one can see that the relative corrections are sensitive to the value of tan (3. 
For tan /3 = 4 and 30, the magnitudes of the corrections can exceed 30% and 20% , 
respectively, for the decay into . For tan/3 = 10, the corrections to the widths of 
the three decay channels are smaller than ones either in the case of tan /? = 4 or in the 
case of tan /3 = 30 . In general, for low tan (3 the top quark contribution is enhanced 
while for high tan j3 the bottom quark contribution become large, and for medium 
tan j3, there are not any the enhanced effects from the Yukawa couplings. So the 
corrections for tan /3 = 4 or 30 are generally larger than those for tan j3 =10, as shown 
in Figs. 3-5 (b). There are some dips and peaks in Figs. 3-5 (b), which arise from the 
singularities at the threshold points m^^ = m^o + rrit and m^^ = wi^^ + mG+{= mw), 
respectively. 

Figs. 6-8 give the tree-level decay widths and the Yukawa corrections as the func- 
tions of niAO for the three decays. We assumed m^-^ = 200GeV, /i = M2 = 200GeV 
and At = Ab = ITeV. The features of the tree level decay widths in Figs. 6-8 (a) are 
similar to Figs. 3-5 (a) , respectively. From Figs. 6-8 (b) we can see that the relative 
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corrections decrease or increase the decay widths depending on tan (3. In most range 
of the mass of ^4°, the relative corrections vary from 27% to 33% for the decay into 
-6% to 20% for the decay into and -9% to -5% for the decay into . There are 
many dips and peaks on the curves in Figs.6-8(b), which come from the singularities 
at the threshold points. For example, at m^o = 235GeV in Fig.8(b) , we have the 
threshold point m^^ = m^o + for tan (3 = 30. 

In Fig. 9 we present the tree- level decay widths and the Yukawa corrections as the 
functions of /i in the case of t2 — ^ + , assuming tan/? = 30, m^^ = 250GeV, 
M2 = lOOGeV, At = 250GeV, Ab = -250GeV and m^o = ISOGeV. In most of the 
parameter fi range , the relative corrections are about from 12% to 32% for the decay 
into h^, and only a few percent for the decay into A^ except near the zero point of 
Fq . For the decay into H^, when /i takes certain values (near about -26 GeV), Fq 
gets very small (< 10~''GeV), and the relative corrections near these values do not 
have a physical meaning. So we cut off the corrections, since perturbation theory 
breaks down here. In order to improve the results, we use the running higgsino mass 
parameter fi{Q) = iJi + 5jl{Q) in the tree-level coupling, and find that the convergence 
of the perturbation expansion becomes better as shown by the dashed line in Fig. 9(b), 
where the region of the parameter /i of breaking down the perturbation theory gets 
smaller (Note that, in fact, the parameter range < 180 GeV has been excluded 



by phenomenology at LEP and Tevatron |]T6|, |3T| ). There are many dips and peaks 
on the curves in Fig.9(b), which come from the singularities at the threshold points. 
For example, at /i = —216 GeV on the solid curve in Fig. 9(b) , we have the threshold 
point = m^o -|- rrit for the decay into H^. 

In Figs. 10-11 we compare the results with the ones presented in an earlier literature 
where the 0{as) SUSY-QCD corrections to the same three decay processes have 
been calculated. We present the tree-level decay widths and the Yukawa corrected 
decay widths as the functions of m^^ and m^o in Figs. 10 and 11, respectively . For 
comparation, we take the same input parameters as in the Ref.p : tan/5 = 3, cos^j = 
0.26, rrii^ = 250GeV, m-g = 600GeV, /i = 550GeV in Figs.10-11, and m^o = ISOGeV 
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in Fig.lO, rrij:^ = 600GeV in Fig. 11. In both Figs, we assumed = 1.12Mq. Our 
numerical results of the tree level decay widths agree with their results except a little 
difference, which is due to the running effects were used in our calculation but not in 
Ref. 1^. The relative corrections in Fig.lO vary from -22% to 26% for the decay into 

-60% to -4% for the decay into H^, and -5% to 0% for the decay into The 
relative corrections in Fig. 11 vary from -1% to 23% for the decay into h^, -24% to 60% 
for the decay into H^, -4% to -1% for the decay into A^. After comparing with Figs. 3 
and 5 in Ref.0, we can see that the Yukawa corrections are comparable to the 0{as) 
SUSY-QCD corrections for the decays into and H^, but smaller than the 0{as) 
SUSY-QCD corrections for the decays into A^. There are two dips at m^o = 348GeV 
and 352GeV on the solid curve of the decay into ho in Fig. 11, which come from the 
singularities at the threshold points m^^ = "^ii + f^no- 

Finally, in Fig. 12 we show the numerical improvement of the Yukawa corrections 
as a function of tan/? in five ways of perturbative expansion: (i) the strict on-shell 
scheme (the dotted line), where the top quark pole mass 175.6GeV, the bottom quark 
pole mass 4.25GeV, the on-shell trilinear coupling At and the higssino mass parameter 
fi were used, (ii) the improved scheme (the solid line), in which the QCD, SUSY- 
QCD, and SUSY-Electroweak running quark masses rhq{Q) and the running trilinear 
coupling At{Q) were used, (iii) the complete improved scheme (the dashed line), in 
which the SUSY-Electroweak running parameter /i was also used as well as the same 
running parameters as in (ii), (iv) the rht{Q) running scheme (the dash-dotted 
line), in which only the running top quark mass was used, and (v) the thbiQ) running 
scheme (the dash-dot-dotted line), in which only the running bottom quark mass 
was used. Here we assumed m^-^ = 250GeV, M2 = 200GeV, At = Ab = 900GeV, 
/i = 200GeV, rriAO = ISOGeV and Mg = 1.5M^ = 1.5M^. One can see that, the 
effect of the running of the top quark mass on the corrections can not be neglected 
for low tan/9(< 10), while the effect of the running of the bottom quark mass is quite 
significant for large tan/?(> 40) . The whole running effects with or without the 
running of the parameter /i both make the convergence of the perturbation expansion 
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much better. The relative corrections approach smoothly -5.0% and 14.3% with the 
increasing tan /3 for the improved scheme and complete improved scheme, as shown 
by the solid line and the dashed line in Fig. 12, respectively. 

In conclusion, we have calculated the Yukawa corrections to the widths of the 
heavier top squark decays into lighter top squarks and neutral Higgs bosons in the 
MSSM. These corrections depend on the masses of the neutral Higgs bosons and 
the lighter or heavier top squark, and the parameter /j.. For favorable parameter 
values, the corrections decrease or increase the tree-level decay widths significantly. 
Especially, for high values of tan (3{—30) or low values of tan/?(=4), the magnitudes of 
the corrections exceed at least 20% for the decay into /i° and which are comparable 
to the 0{as) SUSY-QCD corrections. But for the decay into the corrections are 
smaller and the magnitudes of them are less than 10% in most of the parameter 
space. The numerical calculations also show that using the running quark masses 
and the running trilinear coupling At , which include the QCD, SUSY-QCD, and 
SUSY-Electroweak effects and resume all high order (tan /3)-enhanced effects, can 
vastly improve the convergence of the perturbation expansion. We also discuss the 
effects of the running of the higgsino mass parameter // on the corrections, and find 
that they are significant, too, especially for large tan/3. 
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Appendix A 

The following couplings are given in order 0{ht, hb). 
1. squark - squark - Higgs boson 
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(a) squark - squark - /i° 



-\/2mqh 



q'fq 



up 









Cq, 



for I jJq^jj J ^yP® squarks, respectively. We use the abbreviations 
cos a. a is the mixing angle in the CP even neutral Higgs boson sector, 
(b) squark - squark - 
( 



G% 



-\/2mqhq 



So 



So 



V 



—\/2mqhq 



Co 
So 



So 



(c) squark ~ squark - 

I 



Gl = 

2mw 







(d) squark - squark - 

( 



2mw 







-Aq + 11 ■ 



\ 



Aq- II 



cot [3 
tan /5 



cot/5 1 \ 
tan/9 J 

I 



(e) squark - squark - 

g\ = {G\) 



b _ ir^t\T _ 9 I ml tan P + cot P mt{At cot P + fj,) 



V^mw \ nib(AbtainP + 2mtmb/ sin2P 

(f) squark - squark - G^ 



Gl = {GlY 



— ml mt{At — n cot P) 



V2mw V "^b(ywtan/? - A^) 
2. quark - quark - Higgs boson 







cos/3 1 —ig 



V2 



hb sin P 



9 



^/2 \ sin/3 j ' 2mw 
-nib 



htcosp j ' ^y2mw I "^t 
16 







hn 



' V2 

cos P 



COS P 



Cn 



9 



^2 ^ 1 sin /3 j ' 2mw \ -^b 



ITLt 



^ I sin/3 j ' v/2mH/ [ -'^b 
3. quark - squark - neutralino 



(54) 



N* 
Nt 



(55) 



. J ^''kS 

Here is the 4x4 unitary matrix diagonalizing the neutral gaugino-higgsino mass 
matrix 0, gg. 

4. quark - squark - chargino 



11 — pq Y 



k2 



1,5 _ pg I -^6^/02 
'^ik — ^il 



YtVk2 



(56) 



Here U and V are the 2x2 unitary matrices diagonahzing the charged gaugino 
higgsino mass matrix |0, |30[] . 

5. squark - squark - Higgs boson - Higgs boson 

(a) squark - squark - - Hk (k=l,2) 

( mlSu + mln ^ 



bk 



2V2: 



m 



w 







2mtrnb.\/ 
sin 2/3 



with 



(57) 



Sk = (cos a cos P I sin^ sin a cos P / sin^ /3) 
Tfc = (— sin a sin/3/ cos^ /?, cos a sin/?/ cos^ /?) 
Vfc = (sin(/5 — a), cos(/5 — a)) 



(58) 
(59) 
(60) 



(b) squark - squark - H 

( 



'^55 



hi sin^ P 
hi cos^ P 







"'q 





cos^ /3 
sin2/5 



(61) 
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(c) squark - squark - H - G+ 



56 



w 



I 



\ 



ml tan /? 

2 

t 





cot P 



(d) squark - squark - G - (k=l,2,3) 



^6fe — (^6fe) 



i NT _ 9"^ I m^t SGk + ml TGk 



m 



w 







cot/3 1 
-tan/3 J y 





2mtmb/ sin 2/3 



with 



S'Gfe = (cosQ;/sin/3, sin a/ sin /3 icot/3) 
T'Cfc = (sinQ;/cos/3, — cosq;/cos/3 itan/3) 
^G'ik = (~ cos(/3 — a), sin(/3 — a) — i) 



(e) squark - squark - - H/. (k=l,2,3) 



^kk 



Dlk 





2^^" Dlk 



with 



Dlk — (sin^ a/ sin^ /3, cos^ a/ sin^ /3, cot^ /3) 



2m: 



Z}2, 



2to 





^ml D2u 



with 



D2k = (cos^ aj sin^ /3, sin^ a/ sin^ /3, tan^ /3) 



(f) squark - squark - if" - 



-g^m^ sin2Q; 



D2 



12 



-g^mg sin 2a ^2 



4mf, 



with D2 = -l/cos2/3 
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12 



G 



35 







with Dl = -l/sin^/3 

(g) squark - squark - A" - G° 



/ -g^mg sin 2^ 



4mf, 



D2 



-g^ml sin 2a 



Vmg sin 2/3 



Ami 



(72) 



(73) 



^35 



g m^^ sin 



4mf, 



-g^ml sin 2/3 



4mfi 



(74) 



Finally, we defineG^j = Gf^, and also G^^^-^j^ — 0, when k=l,2,5, i.e. there are no 
A^{G^)h^qq,A^{G^)H^qq, or A^{G^)H+qq couplings. 

Appendix B 



We define q — t and 6, q' the SU (2)^ partner of g, and q" — qlork — 1...4 and q" — q' 
for /c = 5, 6. Then we have 

,2 



5m| 



SZho 



9 



(m2^,m,,m,) 



[ml + m^^ - Ao(m^) - ^(m^) - m^^So - {ml - ml)B^] 



E { o [(^'l - sin^ Owf + sin^ ^H^][2m^ - 2Ao(mJ) - m^i^o] 



-2mleqS\v? 9w{IIl - sin^ dw)Bo}{m%,mg,mq), 



%1 
167r2 



[(m^ cot^ + JTife tan^ /?)(m|^+Gi + -Bi - m^Go) - 2m^m^Go 



{mjj+,mt, mb) + E(<^5)ii(G'*5)iiG'o(rn^+, rnf,, m^.), 



Sqmj cos a , r, ^ o 

^-2m2Go + 5i + Gi) (m^ mi, mt) 



IGTT^m^ sin^ /9 
Sg^mf sin^ a 
167r^m^ cos^ /? 



-2mlGo + Bi + mlGi){ml^,mb, mb) 
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3 

^^Ho = 2 2 -2 A -'^^tGo + Bi+ mHGi){mH nit, rrit) 
IGn^mfy sm p " 

3 - - 

S^f^mJcos^o; , ^ 2^ n 2^ \ / 2 \ 
Sg'^ml sin^ a , ^ ^ 2^ \/ 2 n 

3 



JGH = ~ iQ^2^2 i^i - m^tan/3)(mji?o + Ao{mf,) + mj^+Bi) + m^m^(tan/5 

—3 

■ cot /3)So(m|+, mt, mt) + E(<^5)j«(G'6)«j-Bo(m|+, m^- , m^.) 



+yI^EE^(gI6)..aK), 



5mt 



„2 4 _ -1 



^< = T7r^{EE(Gi)^.(G^r)..i?oK,m,,,mHj - 2<7^E[(I4P + l4r) 

x{m1.Bi + Ao{m'^o) + rn^Bo) + 2mtm^o^Re{al^bll)Bo]{rnf.,rnt, m^c,) 
-2/ E[(l4r + l^nK^i + Mml.) + mlBo) 

k=l 

+2mt'rn^+Re{lli,kll)Bo]{rnl,rnt',rn^+)}, 
^^i^ = T^{EE(4)^i(4'')i^G'o(m|,m.,,m^J + 2/i:[(|aLr + 



20 



x{Bi + rri^^Gi — rn^Go) — 2mtm^o^Re{al^bll)Go]{m'^,,mt,m^c,] 

k=l 

,2 



-2mg,m^+Re{l^i,^kf^)Go]{m-^,mg,,m^+)}, 
SzTT^m^ sm /5 



+ tI^ E E(G'?).^(G'f )..i?o(p^ , m,- ) 

+T7-iEE(G'?2)..^K-j 

OZhOhO = 2 2~' ^^hOHO = 2 2~ 

= Yi2{EE(G'^)l.(Gr),2i?o(p^m,,,,mHJ -2/E[(44; + ^^^^ 

fc=l j fc=l 

-2/ E[(44 + klklDip'B, + Ao(mJ^) + mlBo) 

k=l 

+mj,m-+(44; + lllk{i^)Bo]{p\ mt,, m-+)}, 
50, + 5^1 = -^-^^[j:{^{ml) - S*2(m? )], 



2(m|-m|)^ 
3 2 _ _ 

-■-"^ fc=i 

3 2 _ _ 

ng(p') = - E[4.4,5i(p', rn,, m,-J + fj%Bi{p\ mt, m^J], 
fc=i 

3 2 _ _ _ , 

nj'^d?^) = -^Y.l^bkiilijBo{p'^,mb,miJ +mtli,kljBQ{p^,mt,mf^J], 



167r2 
3 

16^ 



fc=i 

2 



3 - - 

n5'^(p^) = j^Y.i'^bllikljBo{p'^,mb,miJ +mtkl/^^jBo{p'^,mt,miJ]. 
^ fc=i 



Here Aq and i?i are one- and two-point Feynman integrals[^, respectively, and 
Gi = dBi/dp\Go = -dBo/dp'^. 
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Figure 1: Feynman diagrams contributing to supersymmetric electroweak corrections 
to t2 tiHi. Hi, i=l,2,3 correspond to (a) is tree level diagram; (6) — (/) 

are one-loop vertex corrections. In diagram (b) q = t for k = 1...4 and q = b for 
k = 5,6. In diagram (d) and (e) q = t for k = 1,2 and q = b for k = 5,6. In diagram 
if),q = b,t 
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Figure 2: Feynman diagrams contributine to renormalization constants. In diaeram 
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Figure 3: The tree- level decay width (a) of t2 tih^ and its Yukawa corrections (6) 
as functions of mf^ for tan (3 = 4, 10, and 30, respectively, assuming m^o = ISOGeV, 
= Ms = 200GeV, At^A^^ 600GeV and Ma = 1.5Mfj = 1.5M^. 
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Figure 4: The tree-level decay width (a) of ^2 — ^i-f^° and its Yukawa corrections (b) 
as functions of m^^ for tan/^ = 4, 10, and 30, respectively, assuming m^o = 150GeV, 
= M2 = 200GeV, A( = ^6 = eOOGeV and Ma = l.SMj; = 1.5M^. 
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Figure 5: The tree- level decay width (a) of ^2 — and its Yukawa corrections (b) 
as functions of m^^ for tan/^ = 4, 10, and 30, respectively, assuming m^o = 150GeV, 
= M2 = 200GeV, At^Ab^ 600GeV and Ma = 1.5M^ = 1.5M^. 
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Figure 6: The tree-level decay width (a) of ^2 — and its Yukawa corrections (6) 
as functions of m^o for tan/? = 4, 10, and 30, respectively, assuming m^^ = 200GeV, 
= M2 = 200GeV, At^A^^ 600GeV and Ma = l.^M^j = 1.5M^. 
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Figure 7: The tree-level decay width (a) of ^2 — ^i-f^° and its Yukawa corrections (b) 
as functions of ni^n for tan/? = 4, 10, and 30, respectively, assuming rrif^ — 200GeV, 
= M2 = 200GeV, At^Af,^ 600GeV and Ma = l.SM^; = 1.5M^. 
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Figure 8: The tree-level decay width (a) of ^2 and its Yukawa corrections (b) 

as functions of m^o for tan/? = 4, 10, and 30, respectively, assuming m^^ = 200GeV, 
;x = Ms = 200GeV, At ^ ^ 600GeV, and Ma = 1-5M^ = 1.5M^. 
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Figure 9: The tree-level decay width (a) of 4 — and its Yukawa corrections 
(6) as functions of fj,, assuming tan (3 = 30, m^^ = 250GeV, M2 = lOOGeV, At — 
250GeV,Ab = -250GeV, m^o = ISOGeV and Ma = 1-5M^ = 1.5M^. 
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Figure 10: The decay width of ^2 ^ hHf as a function of m^- , assuming tan P — 
3,cos^i = 0.26, rrii^ = 250GeV, = 550GeV, m-g = 600GeV, A = A^, m^o = 150GeV 
and Mf) = 1.12Mq. The sohd hnes correspond to the Yukawa-corrected decay widths, 
The dashed hnes correspond to the tree-level decay widths. 
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Figure 11: The decay width of ^2 — ^ hH^ as a function of m^o, assuming tan/3 = 
3,cos^i = 0.26, rrii^ = 250GeV, fx = 550GcV, m-g = 600GeV, A = A and 
M^, = 1.12Mq. The sohd hnes correspond to the Yukawa-corrected decay widths, 
The dashed hnes correspond to the tree-level decay widths. 
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Figure 12: The Yukawa corrections of t2 — ^ tiH^ as a function of tan (3, assuming 
- 250GeV, M2 = 200GeV, At = = 900GeV, = 200GeV, m^o = 150GeV 
and Mq — 1.5M^ — 1.5M^. The dotted hne corresponds to the corrections using the 
on-shell parameters; the dashed hne corresponds to the corrections using the running 
parameters rht{Q), rhb{Q), At, and fi; the sohd hne corresponds to the corrections 
using the same running parameters except the running jl; the dashcd-dottcd hne to 
the improved resuh only using the running mass rht{Q); and the dash-dot-dotcd hne 
to the improved result only using the running mass rhi,{Q). 
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